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ABSTRACT 

The neutron number spectrum computed by the FASTER code at point detectors in a 
ducted water slab configuration is compared with similar results generated by the 05R 
code. Excellent agreement is found, verifying the quick-running FASTER code for this 
type of shielding analysis. 
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SUMMARY 


The neutron number spectrum computed by the FASTER code at point detectors in a 
ducted water slab configuration is compared with similar results generated by the 05R 
code. Excellent agreement is found, verifying the quick-running FASTER code for this 
type of shielding analysis. 


INTRODUCTION 


One of the major engineering perturbations on the weight of mobile reactor shields 
is additional shielding due to the presence of ducts and voids in the shield configuration. 
Exact calculations of neutron and gamma-ray transport in unducted shields are difficult 
enough and quite time consuming. The effects of ducts in realistic shields is a problem 
which has never been treated thoroughly. 

A method has been tested, for a simplified ducted shield problem, which appears to 
offer promise of significantly reducing computer time. The method is incorporated in 
the computer code FASTER (ref. 1) and has sometimes been referred to as multigroup 
Monte Carlo. It has been used extensively in the analysis of shields for nuclear rocket 
configurations (ref. 2). 

Results of ducted shield calculations with FASTER were compared with results of 
similar calculations with 05R (ref. 3). Excellent agreement was obtained; this constitutes 
a verification of the usefulness of FASTER for the analysis of simply ducted shields. 
FASTER has the advantage of substantial computer time savings over Monte Carlo codes 
such as 05R. 
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ANALYSIS 


The code FASTER was used to compute the differential neutron number flux at four 
point detectors in and near a 10-centimeter-diameter cylindrical duct penetrating water 
slabs of 20- and 100-centimeter thickness, respectively. A point fission neutron source 
was located at the center of one end of the duct. The geometric configuration for the 
20-centimeter-slab calculation is shown in figure 1. 
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In the 100-centimeter-slab thickness calculations the point detectors were located at 
(0, 0, 50), (0, 0, 105), (0, 5, 105), and (0, 10, 105), respectively. 

This ducted slab configuration was used because the geometry and problem definition 
are relatively simple and the 20 - centimeter - slab thickness problem has been solved many 
times at Oak Ridge National Laboratory by means of the 05R neutron Monte Carlo code. 
Oak Ridge solutions were obtained from 1000, 10 000, and 40 000 neutron case histories 
for the 20 -centimeter shield and 2000 histories for the 100-centimeter shield. 

For the 100-centimeter-slab calculations, an importance sampling method known as 
the exponential transformation had to be built into 05R to enhance deep penetration of 
neutrons in the direction of the detectors. Only 2000 case histories were run in 05R for 
the 100 -centimeter slab. 

In the comparisons of the 05R and FASTER results, only the scattered spectra are 
compared (the uncollided contribution is left out), so that the argument of expected agree- 
ment because of the overpowering effect of the large uncollided contributions is not valid - 
the agreement is therefore meaningful. 
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PROBLEM DESCRIPTION 


Geometry 

The problem geometry was very clean, consisting of a water slab (fig. 1) bounded by 
planes at x and y = ±100 centimeters. The slab thickness was assigned the values 20 
and 100 centimeters in the two sets of calculations. The entrance face of the slab was the 
z = 0 plane. The exit faces were the z = 20 and 100 centimeter planes. The water slab 
was penetrated by a right circular cylindrical void duct, with axis normal to the entrance 
and exit faces of the slab, and with a 10 -centimeter diameter. 


Source Definition 

The neutron source was a unit point isotropic source at the origin, emitting 1 neutron 
per second, as is standard in shielding calculations of this type. To conserve histories, 
initial source particle trajectories were defined only over the half-space containing the 
water slab. The neutron energy spectrum was defined by 

n(E) = exp(-E/0. 965) sinh ^2. 29 E 0. 01 < E < 10 MeV 

This spectrum was represented by 21 energy groups. The energy groups spanned the 
range from 0.01 to 10 MeV in 0. 5-MeV increments. 


Detector Locations 

For the 20 -centimeter -thick slab, detectors were located at (0, 0, 10), (0, 0, 21), 
(0, 5, 21), and (0, 10, 21). For the 100-centimeter-thick slab, detectors were located 
at (0, 0, 50), (0, 0, 105), (0, 5, 105), and (0, 10, 105). 


RESULTS AND DISCUSSION 

The neutron number spectra computed with 05R for the 1000- and 40 000-case- 
history, 20-centimeter-slab cases are shown at three of the detectors in figures 2 to 4. 
The corresponding spectra computed with FASTER are also shown in these figures. The 
agreement between the 40 000-case-history 05R results and the FASTER results is excel- 
lent. The decreasing slope of the FASTER spectra at about 2. 3 MeV in energy accurately 
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Figure 2. - Differential flux as function of energy for 
20-centimeter slab at detector L 
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Figure 4. - Differential flux as function of energy for 
20-centimeter slab at detector 4. 
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Figure 3. - Differential flux as function of energy for 
20-centimeter slab at detector 2. 
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Figure 5. - Differential flux as function of energy for 
100-centimeter slab at detector 1. 
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Figure 6. - Differential flux as function of energy for Figure 7, ** Differential flux as function of energy for 

100-centimeter slab at detector 2. 100-centimeter slab at detector 3. 

reflects the window in the oxygen cross section at this energy. 

In carrying out the 05R 100 -centimeter -slab calculations, an importance sampling 
method known as the exponential transformation was built into 05R to enhance deep pene- 
tration of neutrons in the direction of the detectors. Only 2000 case histories were run 
in 05R for the 100-eentimeter slab, but sufficient accuracy was obtained. The FASTER 
results are in good agreement with these 05 R results also. The neutron number spectrum 
comparison is displayed in figures 5 to 7. The agreement is excellent - much better than 
is normally observed in this type of spectral comparison. 

In all the comparisons (figs. 2 to 7) of 05R and FASTER spectra, only the scattered 
spectra are compared, so that the argument of expected agreement due to large uncollided 
contributions is not valid - the agreement of 05R and FASTER is meaningful and consti- 
tutes a verification of FASTER for these ducted shield calculations. 

For these FASTER calculations, the computing time is about 11 minutes per detector 
and was relatively insensitive to slab thickness. The 05R calculations presented in fig- 
ures 2 to 7 required computing times a factor of 2 greater than FASTER for the 1000- 
case-history 20-centimeter-slab calculations, and an order of magnitude greater than 
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FASTER for the 2000-case-history 100-centimeter-slab calculations. The saving in com- 
puter time is the big advantage of FASTER, 


CONCLUSIONS 

Based on the results obtained with FASTER, it appears that the code offers a def- 
inite reduction in computing time in the analysis of ducted reactor shields. If FASTER 
works as well in more practical geometric duct configurations, it may be the first Monte 
Carlo code which can actually be used as a quick-running engineering tool in parametric 
ducted shield analytic studies. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 15, 1969, 

126-15-01-03-22. 
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